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Abstract 

We describe a direct method for monitoring the geometrical dimensions of a syn- 
chrotron beam at the target position for internal target installations. The method 
allows for the observation of the proton beam size as well as the position of the 
beam relative to the target. As a first demonstration of the technique, we present 
results obtained by means of the COSY-11 detection system installed at the cooler 
synchrotron COSY. The influence of the stochastic cooling on the COSY proton 
beam dimensions is also investigated. 
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1 Introduction 



Internal cluster target facilities - such as COSY-11 [1] installed at the cooler 
synchrotron COSY-Jiilich [2] - permit the study of the production of mesons 
in the proton-proton interaction with high luminosity (10 31 cm -2 s" 1 ) in spite 
of very low target densities (~ 10 14 atoms cm -2 ). Such conditions minimize 
changes in the ejectiles' momentum vectors due to secondary scattering in the 
target, and hence facilitate the precise study of reactions with cross section 
values at the nanobarn level. 

An exact extraction of absolute cross sections from the measured data de- 
mands a reliable estimation of the acceptance of the detection system. This 
in turn crucially depends on the accuracy of the determination of the posi- 
tion and dimensions of both the beam and the target. Using the COSY-11 
detection system [1] as illustration, we will describe a method for estimating 
the dimensions of the proton beam, based on the momentum distribution for 
elastically scattered protons, which can be measured simultaneously with the 
investigated reaction. 

At the COSY-11 facility (see Fig. 1) the production of short-lived uncharged 
mesons (rj, u>, r}', or <fi) and hyperons (A and E°) is investigated by means 
of the missing mass technique via the pp — > ppX and pp — > pK + Y reac- 
tions, respectively [3-5]. In this way, the four-momenta of positively charged 
particles are fully determined experimentally. The accuracy of the extracted 
missing mass value determines whether the signal from the given meson is 
visible over a background and depends on the precision of the momentum re- 
construction of the registered protons or kaons. Consequently, it varies with 
the detector resolution, the momentum spread and the geometrical spread of 
the accelerator proton beam interacting with the internal cluster target beam. 
The momentum reconstruction is performed by tracing back trajectories from 
drift chambers through the dipole magnetic field to the target (Fig. 1,[1]), 
which is an infinitely thin vertical line in the ideal case. In reality, however, 
the reactions take place in a region of finite dimensions where the beam and 
the target overlap (see Fig. 2). Assuming an infinitesimal target, our analysis 
shows a smearing out of the momentum vectors and hence of the resolution 
of the missing mass signal. However, if the target and beam dimensions are 
known it is possible to determine the average smearing of the missing mass 
originating from this effect and hence to infer for example the natural width 
of the measured mesons, provided that the beam momentum spread is known 
as well. 
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2 Description of the method 



Part of the COSY-11 detection setup used for the registration of elastically 
scattered protons is shown in Figure 1. Trajectories of protons scattered in the 
forward direction are measured by means of two drift chambers (Dl and D2) 
and a scintillator hodoscope (SI), whereas the recoil protons are registered in 
coincidence with the forward ones using a silicon pad detector arrangement 
(Si) and a scintillation detector (S4). The two-body kinematics gives an un- 
ambiguous relation between the scattering angles 0i and 2 of the recoil and 
forward protons (see Fig. 1). Therefore, as seen in Figure 3, events of elas- 
tically scattered protons can be identified from the correlation line formed 
between the position in the silicon pad detector Si, and the scintillator ho- 
doscope SI, the latter measured by the two drift chamber stacks. For those 
protons which are elastically scattered in the forward direction and deflected 
in the magnetic field of the dipole, the momentum vector at the target point 
can be determined. According to two-body kinematics, momentum compo- 
nents parallel and perpendicular to the beam axis form an ellipse. A section 
of this ellipse is shown as a solid line in Figure 4a superimposed on the data 
which were selected according to the correlation criterion from Figure 3 for 
elastically scattered events. 

In Figures 3 and Figure 4a, the elastically scattered events stand out clearly 
against a low level background. But more importantly the mean of the elasti- 
cally scattered data is significantly shifted from the expected line, indicating 
that the reconstructed momenta are on average larger than expected. This dis- 
crepancy is difficult to explain in terms of the uncertainties of either the proton 
beam momentum or the proton beam angle. For instance, to account for this 
discrepancy, the beam momentum must be changed by more than 120 MeV/c 
(dashed line in Figure 4a), which is 40 times larger than the conservatively 
estimated error of the absolute beam momentum (± 3 MeV/c) [6,7]. Simi- 
larly, the effect could have been corrected by changing the beam angle by 
40 mrad (see dotted line in Figure 4a), which exceeds the uncertainty of the 
beam angle (± 1 mrad [6]) by at least a factor of 40. 

A reasonable agreement is obtained, however, by shifting the assumed reaction 
point relative to the nominal value by —0.2 cm perpendicular to the beam axis 
towards the center of the COSY-ring, along the X-axis defined in Figure 1. The 
experimentally extracted momentum components obtained in this manner are 
shown in Figure 4b, together with the expectation depicted by the solid line. 
There is good agreement and the data are now spread symmetrically around 
the ellipse. This spread is essentially due to the finite extensions of the cluster 
target and the proton beam overlap (about ±0.2 cm as inferred from the shift 
in the target position required). Other contributions, including the spread of 
the beam momentum and still some multiple scattering events, appear to be 
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negligible. We now consider in more detail the influence of the effective target 
dimensions and the spread of the COSY beam. 

We assume the target beam to be described by a cylindrical pipe, with a di- 
ameter of 9 mm, homogeneously filled with protons [8,9]. Furthermore, the 
COSY proton beam density distribution is described by Gaussian functions 
with standard deviations ax and Oy for the horizontal and vertical direc- 
tions (see Fig. 2), respectively [10,11]. Monte-Carlo calculations are performed 
varying the distance between the target and the beam centres (Ax), and the 
horizontal proton beam extension (ax)- 

In order to account for the angular distribution of the elastically scattered pro- 
tons, an appropriate weight w was assigned to each generated event according 
to the differential distributions of the cross sections measured by the EDDA 
collaboration [12]. The generated events were analyzed in the same way as 
the experimental data. A comparison of the data from Figure 4a (data ana- 
lyzed using the nominal interaction point) with the corresponding simulated 
histograms allows both ax and Ax to be determined. For finding an estimate 
of the parameters ax and Ax we construct the x 2 statistic according to the 
method of least squares: 

2 _ {otNf + bj — Nf) 2 
X ^ a 2 J^ w' + Nf + k' {) 

i th bin 



where Nf and Nf — w denote the content of the i th bin of the Pj_-versus- 

i th bin 

P|l spectrum determined from experiment (Fig. 4a) and simulations, respec- 
tively. The background events hi in the i th bin amount to less than one per 
cent of the data [13] and were estimated by linear interpolations between the 
inner and outer part of a broad distribution surrounding the expected ellipse 
originating from elastically scattered protons. The free parameter a allows ad- 
justment of the overall scale of the fitted Monte-Carlo histograms. Thus, by 
varying the a parameter, the Xmin f° r eac h P a i r °f a x and Ax was established 
as a minimum of the x 2 { a ) distribution. 

The results of the Monte-Carlo calculations are shown in Figure 5. 
Figure 5a shows the logarithm of Xmin as a function of the values of ax and Ax- 
The Xmin function has a valley identifying a minimum and allowing the unique 
determination of the varied parameters ax and Ax- The overall minimum of 
the Xmin (°~x, Ax )-distribution, at a x = 0.2 cm and Ax = - 0.2 cm is more 
evident in Figures 5b and 5c, which show the projections of the valley line 
onto the respective axis. 

The same results were obtained when employing the Poisson likelihood x 2 
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derived from the maximum likelihood method [14,15] 



X 2 = 2 • ]T [aN? + h 



Nt + Nt H 



aNf + bi 



)], 



(2) 



The vertical beam extension of ay = 0.51 cm was established directly from 
the distribution of the vertical component of the particle trajectories at the 
centre of the target. This is possible, since for the momentum reconstruction 
only the origin of the track in the horizontal plane is used, while the verti- 
cal component remains a free parameter. As shown in Figure 6, the recon- 
structed distribution of the vertical component of the reaction points indeed 
can be described well by a Gaussian distribution with the a = 0.53 cm. The 
width of this distribution is primarily due to the vertical spread of a proton 
beam. The spread caused by multiple scattering and drift chambers resolution 
was determined to be about 0.13 cm (standard deviation) [13]. Therefore, 



The effect of a possible drift chamber misalignment, i. e. an inexactness of 
the angle f3 in Figure 1, was estimated to cause a shift in the momentum 
plane by 0.15 cm for A X - This gives a rather large systematical bias in the 
estimation of the absolute value of Ax, but it does not affect the parameter ax, 
and still allows for the determination of relative beam shifts Ax during the 
measurement cycle. 



3 Monitoring of the beam geometry during the measurement cycle 

It is possible to use the method described here to control the size and position 
of the beam relative to the target. This was demonstrated during the experi- 
mental run performed in February 1998 [5] , when for the first time longitudinal 
and horizontal stochastic cooling were used at the COSY accelerator. 
The horizontal stochastic cooling [10] is used to squeeze the proton beam in 
the horizontal direction until the beam reaches the equilibrium between the 
cooling and the heating due to the target. The size of an uncooled beam in- 
creases during the cycle, as can be seen in Figure 7a, which shows the spreading 
of the beam in the vertical plane during the 60 minutes cycle. This was ex- 
pected since stochastic cooling in the vertical plane was not used in this case. 
The influence of the applied cooling in the horizontal plane is clearly visible 
in Figure 7b. During the first five minutes of the cycle the horizontal size of 
the beam of about 2 • 10 10 protons was reduced by a factor of 2, reaching 
the equilibrium conditions, and remaining constant for the rest of the COSY 
cycle. 
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Figure 8 depicts the accuracy of the momentum reconstruction, reflected in 
the spread of the data, which is mainly due to the finite horizontal size of the 
beam and target overlap. The upper and lower panel correspond to the first 
and the last minute of the measurement cycle, respectively. The data were 
analyzed correcting for the relative target and beam shifts Ax as determined 
by the method discussed above. The movement of the beam relative to the 
target during the cycle is quantified in Figure 7c. The shift of the beam denotes 
also changes of the average beam momentum, due to the nonzero dispersion at 
the target position. The beam reaches stability after some minutes. This can 
be understood as the equilibration between the energy losses when crossing 
1.6 • 10 6 times per second through the H2 cluster target, and the power of the 
longitudinal stochastic cooling which can not only diminish the spread of the 
beam momentum but also shifts it as a whole. 



4 Conclusions 

We presented a method which allows the proton beam dimensions and its 
position relative to the target to be controlled for internal beam experimen- 
tal facilities. The technique is based on the measurements of the momentum 
distribution of elastically scattered protons and its comparison with the distri- 
butions simulated with different beam and target conditions. We demonstrated 
the applicability of the method in controlling the proton beam during mea- 
surement when stochastic cooling was used at the synchrotron COSY for the 
first time. 
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Fig. 1. Schematic diagram of the COSY-11 detection setup. Only detectors used 
for the measurement of elastically scattered protons are shown. Numbers used on 
the silicon pad detector (Si) and below the scintillator hodoscope (SI) indicate the 
order of segments. Dl and D2 denote drift chambers. The Xgi axis is defined such 
that the first segment of the SI ends at 80 cm and the sixteenth ends at -80 cm. 
The proton beam, depicted by a shaded line, circulates in the ring and crosses each 
time the H2 cluster target installed in front of one of the bending dipole magnets 
of the COSY accelerator. 
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TARGET 




Fig. 2. Schematic depiction of the relative beam and target setting. Seen from above 
(upper part), and from aside (lower part), ax and cry denote the horizontal and 
vertical standard deviation of the assumed Gaussian distribution of the proton beam 
density, respectively. The distance between the centres of the proton and the target 
beam is described as Ax- 
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Fig. 3. Identification of elastically scattered protons from the correlation of hits 
in the silicon detector Si and the SI scintillator hodoscope. Note that the number 
of entries per bin is given in a logarithmic scale, ranging from 1 (smallest box) to 
19000 (largest box). 
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Fig. 4. a) Perpendicular versus parallel momentum components with respect to 
the beam direction of particles registered at a beam momentum of 3.227 GeV/c. The 
number of entries per bin is shown logarithmically. The solid line corresponds to the 
momentum ellipse expected for elastically scattered protons at a beam momentum 
of 3.227 GeV/c, the dashed line refers to a beam momentum of 3.350 GeV/c, and 
the dotted line shows the momentum ellipse obtained for a proton beam inclined 
by 40 mrad. 

b ) The same data as shown in a) but analyzed with the target point shifted by 
-0.2 cm perpendicularly to the beam direction (along the X-axis in Figure 1). The 
solid line shows the momentum ellipse at a beam momentum of 3.227 GeV/c. 
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Fig. 5. a) Xmin as a function of ax and Ax- The number of entries is shown in a 
logarithmic scale, b) Xmin as a function of Ax- c) Xmin as a function of ax- 
The minimum value of Xmin ^ larger than two. However, it decreases when only 
a part of the 60 minutes long COSY cycle is taken into account (see section 3). 
This is due to the fact that the beam changes during the cycle and its form, when 
integrated over the whole cycle time, does not suit perfectly to the form assumed 
in Monte-Carlo simulations. 
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Fig. 6. Distribution of the vertical component of the reaction points determined by 
tracing back trajectories from the drift chambers through the dipole magnetic field 
to the centre of the target (in the horizontal plane). "Tails" are due to secondary 
scattering on the vacuum chamber and were parametrized by a polynomial of second 
order. The solid line shows the simultaneous fit of the Gaussian distribution and 
the polynomial of second order. 
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Fig. 7. The vertical a) and horizontal b) beam width (one standard deviation) 
determined for each five minutes partition of the COSY cycle. 
Note that the beam width ax determined for each five minutes period is smaller 
than the ax averaged over the whole cycle (compare Figure 5c). This is due to the 
beam shift relative to the target during the experimental cycle as shown in the panel 
on the bottom left. 

c) Relative settings of the COSY proton beam and the target centre versus the time 
of the measurement cycle, d) Changes of the luminosity during the cycle. A steep 
decrease on the beginning is caused by the movement of the beam out from the 
target centre. 

The vertical error bars in pictures b) and c) denote the size of the step used in the 
Monte-Carlo simulations (± 0.025 cm ; the bin width of Figures 5b and 5c). 
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Fig. 8. Perpendicular versus parallel momentum components with respect to the 
beam direction of particles registered at a beam momentum of 3.227 GeV/c, as 
measured during the first minute (a), and the last minute (b) of the 60 minutes 
long COSY cycle. The data at the first minute were analyzed with Ax = -0.15 cm, 
and the data of the 60's minute with Ax = -0.25 cm, see figure 7c. The number of 
entries per bin is shown on a logarithmic scale. The solid line corresponds to the 
momentum ellipse expected for protons scattered elastically at a beam momentum 
of 3.227 GeV/c. 
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